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SECTION  1 
INTRODUCTKM 


BACKGROUND 

Since  the  signing  of  the  Limited  Test-Ban  Treaty  in  1963  which 
prohibited  testing  of  nuclear  weapons  in  the  atmosphere,  it  has  been 
necessary  to  test  the  survivability,  vulnerability  and  hardness  of 
military  systems  to  nuclear  weapons  blast  and  shock  effects  through 
simulation  methods  using  conventional  high  explosives.  Since  the  mid- 
1960s,  researchers  have  successfully  adapted  both  military  and 
commercial  explosives  to  the  problems  of  nuclear  weapons  effects  (NWE) 
simulation.  After  25  years  of  testing,  a  wide  variety  of  experiments 
have  been  conducted  using  quantities  of  explosives  ranging  from  a  few 
kilograms  to  several  million  kilograms.  As  the  sophistication  of 
experiments  increased,  the  performance  and  characterization  of  the 
explosives  used  in  these  experiments  became  of  greater  importance, 
expecially  with  the  development  of  large  finite  difference  hydrodynamic 
codes  used  to  design,  predict  and  analyze  major  test  events.  The  trend 
in  recent  years  has  been  toward  larger  scale  tests  in  which  the 
explosive  cost  is  an  important  consideration. 

The  development  of  emulsion  explosives  technology  in  the  past  20 
years  in  the  U.  S.  has  had  a  strong  impact  on  the  commercial  explosives 
industry  and  their  main  customer,  the  mining  industry.  Only  in  recent 
years  has  serious  consideration  been  given  to  military  applications 
of  emulsions.  With  the  wide  variety  of  emulsion  formulations  available, 
the  variety  of  possible  additives  (molecular  explosives,  powdered 
metals,  oxidizers,  etc.)  and  the  control  over  density  (microballoons) 
and  mechanical  properties  (formulation  and  processing),  a  wide  range  of 
detonation  and  mechanical  properties  appear  to  be  potentially  available. 

Of  particular  interest  to  the  NWE  technical  community  is  the  use  of 
emulsions  that  can  be  cast  to  form  rigid  charges  in  a  variety  of  shapes. 
For  instance,  there  is  current  interest  in  the  use  of  very  energetic 
explosives  in  small  spherical  charge  experiments  (500  kg)  and  in  the  use 
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of  Less  energetic  (and  less  expensive)  explosives  in  large  spherical 
charge  experiments  (500,000  kg),  both  for  height-of -burst  experiments. 
The  issues  of  safety,  thermal  stability,  initiation  sensitivity, 
detonation  performance,  mechanical  properties,  processing  methods, 
handling  and  transportation,  environmental  durability  and  cost  are  all 
important . 

Emulsion  explosives  exhibit  certain  inherent  advantages  over  other 
commercial  explosives  such  as  ammonium  nitrate/fuel  oil  (ANFO)  and 
slurries,  both  of  which  have  been  used  extensively  by  the  NWE  community. 
These  advantages  are  a  result  of  the  physical  and  chemical  nature  of  the 
emulsion  structure.  The  salt  solution  (oxidizer)  is  discontinous  and  is 
dispersed  as  small  cells  (less  than  10  /im  diameter)  throughout  the  fuel. 
This  results  in  very  good  intimacy  of  mixing  between  the  fuel  and 
oxidizer,  which  in  turn  results  in  a  matrix  with  a  high  density  and  high 
detonation  velocity  relative  to  ANFO.  In  addition,  the  water-soluble 
salt  solution  is  isolated  from  moisture  by  the  thin  fuel  film  resulting 
in  a  water-resistant  matrix. 

Emulsions  are  classed  as  heterogeneous  with  respect  to  their 
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initiation  behavior  in  response  to  mechanical  shock.  Campbell,  et  al  ’ 
showed  in  a  classic  set  of  experiments  that  homogeneous  explosives  (such 
as  liquid  nitromethane)  and  heterogeneous  explosives  (solid  molecular 
explosives,  slurries,  emulsions,  etc.)  behave  in  fundamentally  different 
ways  during  initiation.  In  homogeneous  explosives  initiation  occurs  due 
to  thermal  effects  caused  by  shock  compression.  After  the  initiating 
shock  has  propagated  a  short  distance  into  the  explosive,  initiation 
occurs  at  the  explosive  surface  where  the  shock  was  initially  applied 
and  a  detonation  front  develops.  The  detonation  front  rapidly  overtakes 
the  shock,  after  which,  steady  detonation  occurs.  In  the  case  of  a 
heterogeneous  explosive,  initiation  occurs  due  to  the  presense  of  hot 
spots  and  the  detonation  fiont  develops  immediately  behind  the  shock. 
Generally,  an  overshoot  in  detonation  velocity  occurs  in  heterogeneous 
explosives  when  initiated  by  strong  shocks;  but,  after  the  detonation 
front  has  propagated  for  some  distance,  the  steady-state  detonation 
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velocity  is  achieved.  Generally,  a  much  more  intense  shock  is  required 
to  initiate  homogeneous  explosives  than  heterogeneous  explosives. 

Due  to  the  inherent  insensitivity  of  most  emulsion  formulations,  it 
is  necessary  to  increase  their  initiation  sensitivity  by  the  use  of 
additives.  The  additives  can  be  either  high  density  inert  particles, 
energetic  particles  (with  an  initiation  sensitivity  greater  than  that  of 
the  emulsion)  or  voids.  In  each  case,  it  is  believed  that  localized 
regions  of  very  high  energy  density  (hot  spots)  are  produced  resulting 
in  high  chemical  reaction  rates  which  lead  to  thermal  initiation.  In 
the  case  of  high  density  particles,  the  increased  energy  density  is  due 
to  shock  interactions  caused  by  the  impedance  mismatch  between  the 
particles  and  emulsion  resulting  in  pressure  induced  hot  spots.  In  the 
case  of  energetic  particles,  it  is  likely  that  shock  interaction  effects 
combine  with  the  effects  of  particle  initiation  to  generate  initiation 
in  the  emulsion.  In  the  case  of  voids,  the  increased  energy  density  is 
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due  to  hydrodynamic  collapse.  Johansson  and  Persson  provide  a  good 
review  of  the  mechanisms  believed  to  be  involved  in  the  hot  spot  model 
of  initiation. 
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Lee  conducted  an  extensive  investigation  of  the  effects  of 

microballoon  size  and  loading  on  the  detonation  properties  of  emulsion 
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explosives.  His  research  was  preceded  by  that  of  Yoshlda,  et  al  ,  and 
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Hattorl,  et  al  .  The  main  results  of  Lee's  work  are  siimmarlzed  as 
follows:  (1)  At  a  given  diameter,  as  the  emulslon/mlcroballoon  density 
Increases ,  the  detonation  velocity  also  Increases  until  a  maximum  Is 
reached,  after  which,  the  velocity  begins  to  decrease  until  detonation 
failure  occurs;  (2)  Failure  diameter  Increases  as  the  density  Increases; 
(3)  At  a  given  diameter,  detonation  velocity  Increases  as  the  diameter 
of  the  microballoons  decreases;  and  (4)  Sensitivity  to  shock  Initiation 
Increases  as  the  density  decreases. 

Lee's  results  for  unsleved  microballoons  are  summarized  In  a 
qualitative  manner  In  Figure  1  where  detonation  velocity,  normalized  to 
the  Infinite  diameter  detonation  velocity  of  the  pure  emulsion  (0/D^), 
Is  plotted  as  a  function  of  density,  normalized  to  the  pure  emulsion 
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NORMALIZED  DETONATION  VELOCITY 


NORMALIZED  DENSITY 


Figure  1.  Detonation  velocity  as  a  function  of  density  for  emulsion/ 
microballoon  and  molecular/emulsion  explosive  systems,  normalized 
parameters . 


4 


density  (p/po)-  Thus,  the  region  to  the  left  of  p/Po“l  is  for 
emulsion/microballoon  mixtures.  Thf;  effects  of  explosive  density  and 
diameter  on  detonation  velocity  and  :.ailure  diameter  are  revealed.  In 
reality,  the  detonation  velocity  at  failure  for  a  given  diameter  is 
difficult  to  define  experimentally  and  large  uncertainties  exist.  The 
failure  boundary  in  Figure  1  is  thus  to  be  taken  as  highly  idealized, 
but  nevertheless ,  descriptive  of  the  qualitative  behavior.  The  results 
of  Becker -Kistiakowsky-Wilson  (BKW)  calculations  are  also  shovm  to 
provide  an  idealized  prediction  of  detonation  velocity  for  an  infinite 
diameter  charge.  Note  that  the  behavior  shown  in  Figure  1  is 
descriptive  of  (1)  and  (2)  in  the  previous  paragraph. 

If  a  particulate  molecular  explosive  was  added  to  the  pure  emulsion 
then  an  extension  of  the  results  in  Figure  1  could  be  made  into  the 
region  where  p/pq>\,  assuming  that  the  density  of  the  molecular 
explosive  is  greater  than  that  of  the  pure  emulsion.  In  this  case,  the 
detonation  velocity  at  infinite  diameter  is  expected  to  increase  as  the 
density  of  the  mixture  Increases  with  the  maximum  obtained  when  the 
mixture  reaches  the  limit  of  100  percent  molecular  explosive.  If  the 
failure  diameter  and  velocity  at  failure  of  the  molecular  explosive  is 
known,  then  it  is  possible  to  construct  a  qualitative  prediction  of  the 
effects  of  molecular/emulsion  explosive  density  and  charge  diameter  on 
detonation  velocity  and  failure  diameter.  This  qualitative  prediction 
is  shown  in  Figure  1  where  the  nominal  values  of  octahydro-1,3,5,7- 
tetranitro-l,3,5,7-tetrazocine  (HMX)  density  and  detonation  velocity 
have  been  assumed  for  the  molecular  explosive.  Note  that  at  a  given 
charge  diameter,  the  path  of  constant  charge  diameter  through  D/D^^  vs 
p/Po  space  is  continous  if  mixtures  of  all  three  components 
(microballoons,  emulsion,  molecular)  are  not  considered  and  only 
diameters  greater  than  the  failure  diameter  of  the  pure  emulsion  are 
allowed.  If  diameters  less  than  that  of  the  failure  diameter  of  the 
pure  emulsion  are  allowed,  then  it  is  observed  that  the  lines  of 
constant  diameter  are  dlscontinous,  being  "cut  off"  by  the  detonation 
failure  boundaries.  Note  again  that  the  purpose  of  Figure  1  is  to 
represent  only  the  qualitative  behavior  of  the  mlcroballoon/emulsion  and 
molecular/emulsion  mixtures. 
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It  would  be  a  monumental  task  to  conduct  an  experimental  program  in 
which  the  entire  three -component  (microballoon/emulsion/molecular) 
system  was  characterized.  In  addition,  it  is  not  clear  that  there  would 
be  applications  where  the  detonation  characteristics  resulting  from  a 
particular  mixture  of  all  three  components  would  be  of  interest. 
Regardless,  the  portions  of  the  system  represented  in  Figure  1 
(microballoon/emulsion  and  molecular/emulsion  mixtures)  are  of  obvious 
interest  and  a  complete  characterization  is  warranted.  Thus,  it  is  the 
objective  of  this  research  program  to  investigate  the  detonation 
characteristics  of  a  molecular/emulsion  explosive  system. 

OBJECTIVE 

The  objective  of  this  research  is  to  characterize  the  detonation 
performance  of  molecular /emulsion  explosive  mixtures.  Specifically, 
this  research  will  generate  the  following  experimental  data;  (1) 
detonation  velocity  as  a  function  of  molecular  explosive  loading  and 
charge  diameter;  and  (2)  failure  diameter  and  detonation  velocity  at 
failure  as  a  function  of  molecular  explosive  loading.  These  results 
will  be  used  to  verify  the  qualitative  predictions  shown  on  the  right 
side  of  Figure  1. 

APPROACH 

HMX  was  selected  as  the  molecular  explosive  to  be  used  in  this 
research  program.  This  selection  was  made  for  two  reasons.  First,  a 
very  energetic,  well -characterized  molecule  was  required.  HMX  clearly 
satisfies  this  requirement.  Second,  the  Air  Force  Weapons  Laboratory 
(AFWL)  had  HMX  available  to  support  this  program.  Since  the  time  for 
this  research  was  restricted,  material  availability  was  an  Important 
consideration. 

The  emulsion  used  in  this  program  was  prepared  at  the  Center  for 
Explosive  Technology  Research  (CETR) ,  New  Mexico  Institute  of  Mining  and 
Technology  (NM  Tech);  however,  HMX/emulsion  mixing  was  performed  at  the 
Los  Alamos  National  Laboratory  (LANL) .  This  was  necessary  because  the 
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materials  processing  laboratory  construction  at  CETR  had  not  progressed 
to  the  point  that  HMX  drying  and  mixing  could  be  performed.  All 
experimental  work  was  performed  at  the  CETR  Eagle  Site  located  west  of 
the  NM  Tech  campus . 
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SECTION  2 

EXPERIMENT  DESCRIFTIW 


The  primary  data  collected  were  from  cylindrical  and  conical  charge 
experiments.  Charges  were  prepared  by  filling  high-dersity  urethane 
foam  charge  containers  with  the  HMX/emulsion  mixes.  Detonation  front 
position- time  data  were  obtained  with  ionization  pins,  shorting  pins  and 
streak  camera.  In  this  section,  test  articles,  HMX/emulsion  preparation 
procedures,  test  article  loading  procedures  and  instrumentation  systems 
are  described. 

CHARGE  DESIGNS 

Cylindrical  Charges.  Rigid  high-density  :.losed-cell  urethane  foam 
provided  by  General  Plastics  Inc.,  Tacoma,  VA,  was  selected  for  use  in 
fabrication  of  the  cylindrical  charge  containers.  This  material  is 
available  at  several  densities.  A  density  of  0.24  g/ml  was  chosen  as  a 
compromise  between  minimum  cell  size  (high  density)  and  minimum 
confinement  (low  density).  A  small  cell  size  is  desirable  to  provide  a 
smooth  machined  surface.  This  material  is  machined  easily  with  either 
wood  or  metal  cutting  techniques. 

The  cylindrical  charge  containers  were  designed  to  provide  a  run-up 
length  of  10  diameters  prior  to  the  first  time -of- arrival  (TOA)  pin. 
Four  TOA  pins  were  used  in  each  charge  with  a  spacing  between  pins  equal 
to  one  charge  diameter  or  20  mm,  whichever  was  greatest.  The  minimum 
spacing  of  20  mm  was  imposed  to  Insure  that  the  combined  effects  of 
uncertainties  in  pin  spacing  (0.05  mm)  and  instrumentation  sample  rate 
(100  MHz)  would  always  provide  less  than  0.5%  uncertainty  in  the 
detonation  velocity  calculated  from  the  TOA  data.  The  pins  were 
oriented  normal  to  the  charge  axis  and  protruded  a  maximum  of 
approximately  1  mm  into  the  charge  to  minimize  perturbations  to  the 
detonation  front.  Even  though  some  curvature  to  the  detonation  front  is 
known  to  be  present,  for  steady-state  detonation  the  velocities  at  the 
center  of  the  charge  and  at  its  outer  radius  are  identical. 
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For  Che  limited  nature  of  this  research  program,  cylindrical  charge 
containers  for  only  five  diameters  were  fabricated.  Table  1  and 
Figure  2  summarize  the  important  design  features  of  the  cylindrical 
charge  containers.  Each  container  was  fabricated  in  two  segments,  as 
indicated  in  Figure  2,  to  facilitate  machining  and  to  maintain 
dimensional  tolerances  on  the  bore.  For  the  two  smallest  charges  the 
surfaces  of  the  bores  were  sealed  with  paraffin  to  insure  a  smooth 
explosive -urethane  interface.  This  was  felt  to  be  necessary  for  the 
smaller  bores  where  the  urethane  cell  size  was  on  the  order  of  2%  of  the 
bore  diameter. 

As  shown  in  Figure  2(b),  small  1 . 6-mm-diameCer  acrylic  rods  were 
placed  in  the  sides  of  Che  3.18-  and  6 . 35-mm-diamecer  charge  containers 
Co  provide  for  light  transmission  in  streak  camera  experiments.  This 
technique  was  used  for  the  smaller  charges  to  avoid  perturbation  of  the 
detonation  front  by  the  protrusion  of  the  ionization  pins  into  the  bore. 
The  acrylic  rods  were  positioned  flush  with  the  bore. 

Additional  cylindrical  charges  were  prepared  for  plate  dent 
experiments.  These  consisted  of  cardboard  and  PVC  pipe  containers 
positioned  on  the  center  of  steel  plates.  No  TOA  measurements  were  made 
in  the  plate  dent  experiments . 

Conical  Charges.  The  conical  charges  were  fabricated  through  a 
cast  operation  using  0.08  g/ml  urethane  foam.  General  Plastics  Inc. 
fabricated  the  molds  and  performed  the  casting  operations.  Four 
sections  of  different  sizes  were  fabricated.  Detailed  descriptions  of 
the  sections  and  their  dimensions  are  provided  in  Figure  3  and  Table  1. 
The  different  size  sections  permitted  selection  of  combinations  which 
would  bracket  the  expected  failure  diameter  for  any  given  HMX/emulsion 
mixture.  Because  the  urethane  forms  a  continous  skin  at  the  mold 
boundaries,  the  cast  procedure  resulted  in  a  conical  charge  container 
with  a  mirror -smooth,  precision  surface.  The  sections  were  assembled  by 
aligning  them  with  a  short  tapered  plug  sized  to  fit  the  mating  ends. 
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(a)  Nominal  Values.  See  Figures  2  and  3  for  nomenclature. 


I 
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(b)  STREAK  CAHERA  CONFIGURATION 
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Figure  3.  Conical  charge  container. 


Note  in  Figure  3  that  a  single  segment  of  a  cylindrical  charge 
container  is  attached  to  the  large  end  of  the  conical  charge  container 
to  serve  as  a  driver  stage  to  insure  that  the  detonation  front  is  near 
steady- state  conditions  as  it  enters  the  conical  section.  Shock 

position- time  data  was  measured  using  a  resistance  probe  positioned 
along  the  centerline  of  the  charge  for  the  section  I -II  conical  charges 
(diameter  from  50.8  to  12.7  mm).  The  two  smaller  conical  charges 

assembled  from  sections  II-III  and  III-IV  were  configured  for  streak 
C5imera  measurements  to  avoid  what  was  felt  would  be  excessive 
perturbation  by  the  resistance  probe.  This  was  accomplished  by 
machining  a  narrow  slot  along  the  length  of  the  conical  sections  and 
inserting  an  acrylic  strip  flush  with  the  conical  (interior)  surface. 
The  acrylic  strip  served  to  provide  immediate  quenching  of  the 
detonation  light  to  provide  a  well-defined  image  of  the  detonation  front 
trajectory.  The  camera  viewing  slot  was  1.6  mm  wide  and  aligned  along 
the  center  of  the  acrylic  strip. 

HHX/QftTLSKMI  PREPARATION 

Emulsion  Preparation.  The  emulsion  formulation  selected  was  as 
follows : 


Salt  Solution  Ammonium  Nitrate  (AN)  71.5% 

Sodium  Nitrate  (SN)  10.0% 

Distilled  Water  12.0% 


Fuel 


Light  Mineral  Oil/Emulsifier  6.5% 


The  salt  solution 
cells  of  from  1-  to 
phase) .  The  fuel 
emulsifier. 


(discontinous  phase)  is  dispersed  as  tiny  spherical 
lO-pa  diameter  throughout  the  fuel  (continuous 
consists  of  2.7  parts  light  mineral  oil  to  one  part 


The  fuel  was  prepared  by  adding  the  emulsifier  components  to  the 
light  mineral  oil  at  a  temperature  of  80  and  stirring  for  4  h.  The 
emulsifier  components  consisted  of  100  parts  of  polyisobutenyl  succinic 
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anhydride  (PSA)  and  5.62  parts  of  a  50/50  mixture  of  diethanolamine  and 
monoethanolamine .  This  formulation  was  developed  based  on  discussions 

7 

with  Griffith  and  Oxley  . 

The  emulsion  was  prepared  as  follows.  The  salt  solution 
components  were  mixed  together  cold  and  then  heated  to  107  °C.  Complete 
melting  (as  indicated  by  a  clear  solution)  occurred  at  82  °C.  The  fuel 
was  heated  in  an  oil  bath  to  107  °C  and  then  placed  in  a  preheated 
stainless  steel  thermally  insulated  mixing  bowl  (Kitchen  Aid) .  The  salt 
solution  was  added  to  the  fuel  (all  at  once)  and  mixed  at  a  moderate 
speed  (#4)  for  3  min  by  a  Kitchen  Aid  blender  using  a  wire  whip.  The 
mixing  speed  was  then  increased  to  #10  for  3  min.  Typically,  after 
approximately  1  min  at  the  high  speed,  the  mixture  would  suddenly 
emulsify.  The  emulsion  was  then  placed  in  a  plastic  mixing  bowl  and 
beat  at  high  speed  with  a  kitchen  mixer  for  an  additional  3  min, 
stopping  every  45  seconds  to  scrape  the  sides  of  the  bowl.  Temperature 
at  completion  of  mixing  was  typically  88  °C,  well  above  the 
crystalization  temperature  of  the  salt  solution  (82  °C>. 

Numerous  samples  were  taken  for  density  measurements  and 
examination  under  an  optical  microscope.  Due  to  air  entrainment  during 
mixing,  a  good  measure  of  the  void-free  emulsion  density  was  not 
obtained.  It  was  not  possible  to  remove  the  entrained  air  by  vacuum 
due  to  the  stiffness  of  the  emulsion.  Techniques  were  developed  later, 
after  HMX/emulslon  mixing,  to  remove  most  of  the  voids.  A  good  measure 
of  the  void- free  emulsion  density  was  finally  obtained  by  adding  a  known 
quantity  of  light  mineral  oil  to  reduce  the  stiffness  of  the  emulsion 
which  allowed  the  entrained  air  to  escape.  By  measuring  the  density  of 
the  diluted  emulsion  and  considering  the  density  of  the  added  mineral 
oil,  the  void-free  density  was  determined  to  be  1.43  g/ml. 

Figure  4  shows  a  typical  microscope  photograph  of  the  emulsion  at  a 
power  of  160X.  Photographs  were  also  taken  at  160X  of  10/90 
HMX/emulslon  mixtures  to  determine  the  effects  of  additional  mixing  of 
the  emulsion  matrix  at  room  temperatures  when  loaded  with  a  particulate. 
The  emulsion  cell  size  distribution  was  measured  for  each  case  and  the 
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results  in  terms  of  frequency  and  mass  were  plotted  in  Figure  5.  The 
frequency  and  mass  are  normalized  but  the  results  can  be  directly 
related  to  the  actual  frequency  and  mass  distributions  over  the  cell 
sizes.  The  results  in  Figure  5  are  as  expected.  The  nominal  mass  mean 
cell  size  of  the  pure  emulsion  is  about  7  /xm,  and  the  effects  of  mixing 
in  the  HMX  is  to  reduce  the  mass  mean  cell  size  of  the  emulsion  to 
around  4  /im.  It  is  assumed  that  higher  loadings  of  HMX  in  the  emulsion 
would  result  in  even  smaller  cell  sizes. 

HMX/Enulslon  Mixing  Operations.  The  original  plan  for  HMX/emulsion 
mixing  was  to  produce  mixtures  ranging  from  10%  to  70%  HMX  using  a  blend 
of  HMX  grades  composed  of  35%  fine  and  65%  coarse  grade  HMX.  This 
particular  blend  was  selected  because  (1)  the  grades  necessary  for  this 
blend  were  available  and  (2)  this  was  the  blend  selected  by  the  Defense 
Nuclear  Agency  in  the  casting  of  500-kg  HMX  charges  using  a  propellant 
binder  (85%  HMX/15%  binder).  A  50/50  mix  using  only  coarse  grade  HMX 
(50C/50)  and  another  using  only  fine  grade  KMX  (50F/50)  were  planned  to 
investigate,  in  a  veiry  preliminary  way,  the  effects  of  particulate  size 
on  detonation  performance.  As  will  be  shown,  this  plan  was  revised 
immediately  after  mixing  operations  began. 

The  HMX  blend  was  used  initially  in  a  40/60  mix  resulting  in  a  very 
stiff  matrix  with  the  consistency  of  a  soft,  nontacky  putty.  Next,  the 
50F/50  mix  was  attempted,  and  the  matrix  apparently  failed,  resulting  in 
a  dry,  lumpy  noncohesive  mix.  A  l-kg  batch  of  50C/50  was  prepared  in  a 
small  blender  resulting  in  a  mix  with  the  consistency  of  a  stiff  paste. 
It  was  then  concluded  that  the  fine  grade  HMX  was  presenting  the 
emulsion  (in  particular,  the  continous  phase  of  the  emulsion)  with 
excessive  surface  area.  The  decision  was  made  to  use  only  coarse  grade 
HMX  for  the  remaining  mixes .  The  terminology  was  changed  so  that  the 
coarse  grade  HMX  had  no  identifying  letter  (such  as  "C")  and  the  blend 
was  identified  by  the  letter  B  in  the  specification  of  the  HMX/emulsion 
mixture  ratio,  i.e.,  40B/60. 

The  remaining  mixes  were  prepared  using  only  coarse  grade  HMX.  The 
30/70,  20/80  and  10/90  mixes  provided  smooth  matrices  with  decreasing 
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NORMALIZED  FREQUENCY,  MASS  FACTOR 


CELL  DIAMETER  « MICROMETERS) 


Figure  5.  Eiaulslon  cell  size  mass  and  frequency  distribution  for  100% 
emulsion  and  10/90  HMX/emulsion  mixture. 


stiffness.  The  70/30  mix  had  the  consistency  of  a  slightly  damp 
noncohesive  powder.  A  large  quantity  of  50/50  mix  (several  kg)  was 
prepared  resulting  in  the  stiff  paste  as  observed  previously  with  the 
1-kg  50/50  mix.  However,  within  two  days  this  mix  solidified.  It  is 
suspected  that,  due  to  the  high  HMX  loading,  excessive  shear  during 
mixing  resulted  in  matrix  failure  and  crystalization  of  the  salt 
solution.  A  60/40  mix  was  not  attempted  due  to  a  shortage  of  emulsion 
and  the  strong  probability  that  a  successful  mix  would  not  result. 

A  second  mixing  operation  was  conducted  to  (1)  provide  the  40/60 
mix  not  obtained,  during  the  Initial  mixing  operation;  (2)  provide 
additional  quantities  of  20/80  and  50/50  mixes,  and  (3)  attempt  to 
obtain  a  good  60/40  mix.  Care  was  taken  this  time  to  apply  minimum 
shear,  especially  to  the  50/50  and  60/40  mixes.  A  good  40/60  m'xture 
was  obtained,  however,  the  50/50  mix  solidified  as  before  and  the  60/40 
mix  produced  a  dry  but  somewhat  cohesive  consistency,  l.e.,  it  could  be 
packed  to  form  a  relatively  stable  mass.  Thus,  the  experimental  program 
proceeded  with  an  insufficient  quantity  of  50/50  mix  to  conduct  a  51.5- 
mm-dlameter  cylindrical  charge  experiment  planned  as  part  of  the  test 
matrix. 

Coarse  grade  HMX  size  specifications  are  shown  in  Table  2.  Table  3 
presents  the  results  of  the  hazards  analysis  testing  performed  by  LANL 
for  HMX/emulslon  mixtures.  A  notable  result  of  the  data  in  Table  3  is 
the  greatly  reduced  spark  sensitivity  for  HMX/emulsion  mixtures  relative 
to  pure  HMX.  HMX/emulsion  mixtures  of  greater  than  50%  HMX  show 
increased  sensitivity  to  impact,  though  somewhat  less  so  than  pure  HMX. 

HHX/Emulslon  Density  Kessureaents .  Density  measurements  were  made 
by  loading  a  small  cup  with  the  HMX/emulsion  mixtures  using  the  same 
techniques  that  were  used  in  loading  the  charge  containers.  Details  of 
the  loading  procedures  will  be  discussed  under  CHARGE  LOADIHG.  The 
results  are  shown  in  Figure  6  where  density  is  plotted  as  a  function  of 
percent  HMX.  The  data  plotted  for  the  40%  HMX  mixture  are  for  the  40/60 
mix.  The  density  of  the  40B/60  mix  was  essentially  the  same  as  that  of 
the  40/60  mix.  The  density  for  the  60/40  mix  lies  slightly  above  the 
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Table  2.  HMX  specif ications . 


Cltts  2,  Grade  «,  not  araatar  than  2X  RDX 


SIEVE  SIZE  (un) 

PERCENT  PASSING 

297 

100 

125 

98 

44 

75 

Table  3.  Hazards  analysis  data. 


FOKNULAT ION 


SENSITIVITY 

75/25 

50/50 

25/75 

0/100 

lapact  Sanaltivity 

Type  13/13* 

35.6/100 

43.5/95.5 

320/>320 

>320/>320 

OTA 

Exothar*  Peak  <C) 

236 

236 

240 

240 

Onaat  of  Exothara 

<C) 

200 

210 

200 

250 

Pyrolysla  <C) 

256 

258 

264 

270 

Spark  Sanaltivity 

3  aft/18  ail  (J) 

1.73/4.1 

— 

— 

1  .90/8.4 

Nankin  Critical  Taap 

(C) 

209 

•  •  • 

•  •  • 

300 
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curve,  probably  because  it  had  to  be  packed  into  the  charge  container 
which  apparently  drove  some  of  the  entrained  air  out.  The  density  of 
the  70/30  mix  lies  well  below  the  other  data  because  the  noncohesive 
nature  of  that  mix  made  it  impossible  (at  reasonable  pressures)  to  pack 
the  mix  to  form  a  mechanically  stable  mass. 

The  mixture  equation  which  calculates  the 
mixture  of  n  components  of  densities  pj,  Pi>  ■■ 

Xx,  Xi,-  X^  Is 

.  _  _ 1 _ 

^m  Xi/Pi  +  X2/P2  * 

The  curve  shown  in  Figure  6  is  a  plot  of  equation  (1)  where  the 
emulsion  density,  p^,  has  been  adjusted  to  provide  the  best  fit  to  the 
data.  The  emulsion  density  which  gave  the  best  fit  was  pj-  1.389  g/ml. 
The  standard  density  of  1.89  g/ml  for  HMX,  pj,  was  used.  Recall  that  an 
emulsion  density  measurement  was  made  where  all  of  the  entrained  air  was 
removed  to  obtain  a  density  of  1.43  g/ml.  This  measurement,  when 
compared  to  the  value  determined  by  the  fit  of  Equation  (1)  to  the  data, 
Implies  that  the  HMX/emulslon  mixing  process  consistently  Introduced  air 
Into  the  emulsion  at  a  constant  volume  fraction  of  approximately  3%. 

CHARGE  LOADING 

Various  techniques  were  developed  for  loading  of  the  test  articles 
because  of  the  range  of  matrix  stiffnesses  encountered  for  the  various 
HMX/emulslon  mixes.  The  smaller  charge  containers  were  loaded 
successfully  with  10/90,  20/80,  30/70  and  40/60  mixtures  using  manually 
pressed  frosting  bags  with  a  simple  conical  tip.  For  loading  of  the 
larger  charges,  a  simple  pressure  vessel  was  constructed  which  forced 
the  mixture  Into  the  containers  under  air  pressure  produced  by  a  hand 
operated  bicycle  pump.  These  loading  procedures  worked  quite  well  and, 
with  only  one  notable  exception,  resulted  In  uniform  void-free  charges. 


density,  p^,  of  a 
.  p^  and  mass  fractions 
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Due  to  the  very  stiff  nature  of  the  40B/60,  50/50  and  60/40 
mixtures,  a  completely  different  approach  to  test  article  loading  was 
needed.  A  procedure  which  worked  very  well,  even  though  somewhat 
tedious,  was  to  drop  small  balls  of  the  mixture  approximately  one-half 
the  diameter  of  the  bore  and  tamp  it  with  a  wooden  dowel  rod.  For  the 
conical  charges,  a  similar  procedure  was  used  by  loading  from  the  large 
end  using  small  balls  and  dowel  rods  initially  and  increasing  their 
diameters  as  the  conical  container  was  filled. 

Prior  to  filling  a  frosting  bag  or  the  pressure  vessel,  the  mixture 
was  smoothly  spread  over  the  surface  of  a  one-quarter  inch  thick  acrylic 
sheet  using  a  rubber  spatula.  The  mixture  was  worked  to  a  thickness  of 
approximately  1  to  2  mm.  Visible  clumps  of  HMX  were  removed  (>0.5  mm), 
and  then  the  mixture  was  recovered  from  the  acrylic  sheet  with  a  rubber 
spatula  using  a  slow  snowplow  action.  This  procedure  removed  most  of 
the  larger  air  bubbles  entrained  in  the  mixture . 


INSTRDMEinATlOif  SYSTEMS 

Electronic  Measurements.  The  Big  Eagle  Facility  is  configured  for 
electronic  measurements  of  TOA  from  ionization  and  shorting  pins,  and 
voltage  versus  time  measurements  provided  by  resistance  probes.  A 

LeCroy  Model  6900-1  Digital  Oscilloscope  and  Instrumentation  Control 

s 

System  is  the  heart  of  the  measurement  system.  An  IBM  PC  with  CATALYST 
software  is  used  to  set  up  recording  channels  and  analyze  data.  For  TOA 
measurements,  a  pulse-forming  network  provides  the  signal  conditioning 
for  recording  on  the  LeCroy  system.  For  resistance  probe  measurements, 
a  constant  current  power  supply  is  used  to  drive  the  resistance  probe. 
Permanent  digitized  data  records  are  provided  on  disks. 

Screak  Camera.  Streak  camera  data  were  recorded  on  a  Beckman  & 

9 

Whitley  Model  770  rotating  mirror  c«uneca.  The  camera  records  during 
1/12  of  the  revolution  of  the  mirror  on  a  200-mm  length  of  70-mm  film. 
Rotor  position  and  velocity  are  monitored  by  means  of  a  pickup  coil  and 
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armature  placed  over  a  magnetized  rotor  shaft.  The  rotor  is  driven  by  a 
dry  nitrogen  turbine  system. 

Synchronization  and  control  is  accomplished  by  operator  controls 
and  synchronization  electronics.  A  trigger  pulse  is  produced  which  will 
initiate  the  event  at  some  specified  time  increment  in  advance  of  the 
time  at  which  the  mirror  image  encounters  the  film  leading  edge.  The 
system  automatically  produces  the  trigger  pulse  when  the  firing  button 
is  depressed  and  the  turbine  is  up  to  speed. 
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SECTION  3 

THEORETICAL  CALCULATKMS 


Theoretical  calculations  were  performed  for  HMX/emulsion  mixtures 

over  the  range  of  from  100%  emulsion  to  100%  HMX.  Calculations  were 

performed  using  the  TIGER  code  with  the  BKWR  equation  of  state  which 

1  0 

utilizes  the  revised  BKW  parameters  developed  by  Mader  in  fitting 
calculational  results  to  experimental  data  for  RDX.  Since  the 

thermodynamic  modeling  in  TIGER  gives  no  consideration  to  geometry  and 
size,  the  results  are  essentially  representative  of  those  for  an 
infinite  diameter  charge  for  planar  detonation.  Thus,  comparisons  of 
theoretical  and  experimentally  determined  detonation  velocities  must  be 
made  based  on  extrapolations  of  experimental  results  to  infinite 
diameter . 

Mixture  densities  used  In  the  calculations  are  based  on  the 
experimental  results  shown  in  Figure  6.  The  TIGER  library  did  not 
contain  constituent  data  (molecular  composition,  enthalpy  of  formation, 
molar  volume,  state  reference,  specific  heat)  for  light  mineral  oil  and 
the  emulsifier,  PSA,  so  substitutions  were  made  from  the  existing 
library.  Paraffin  is  an  excellent  representation  of  light  mineral  oil 
and  data  for  a  standard  commercially  available  emulsifier.  Span  80,  was 
used  as  an  approximation  for  PSA. 

The  constituents  for  the  calculations  are  summarized  in  Table  4. 
The  results  of  the  calculations  are  summarized  in  Table  5. 

The  oxygen  balance  calculations  in  Table  5  were  performed 
separately  from  the  TIGER  calculations  and  assume  that  all  available 
oxygen  reacts  to  form  CO,  and  H^O.  A  small  code  written  in  BASIC  called 
BALANCE  was  developed  to  perform  the  calculations.  A  listing  of  BALANCE 
is  provided  in  the  appendix. 

The  calculational  and  experimental  results  will  be  compared  in 
Section  4. 
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Table  4.  HMX/emulslon  constituents  for  TIGER  calculations. 
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Table  5.  Results  of  HMX/emulsion  TIGER  calculations. 
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SECTION  4 

RESULTS  AMD  DISCUSSION 

SUMMARY  OF  EXPERIMENTAL  RESULTS 

The  experimental  results  are  summarized  in  Tables  6  through  9.  In 
Table  6  the  results  are  shown  in  terms  of  whether  or  not  successful 
detonation  was  achieved  for  each  cylindrical  charge  experiment.  Table  7 
presents  the  same  results  for  the  conical  charge  experiments . 

Table  8  summarizes  the  experiments  conducted  at  the  Big  Eagle 
Facility  and  Table  9  summarizes  those  conducted  at  the  Little  Eagle 
Facility.  Experiments  designed  for  electronic  TOA  and  resistance  probe 
measurements  were  conducted  at  Big  Eagle.  Experiments  designed  for 
streak  camera  photography  were  conducted  at  Little  Eagle. 

Most  of  the  experiments  were  straightforward  and  the  results  are 
adequately  described  in  the  tables.  For  a  few  experiments,  additional 
remarks  are  required  as  follow: 

BE87-118.  This  was  a  conical  charge  experiment  for  the  20/80 
mix  in  which  resistance  probe  measurements  were  made  to 
determine  detonation  velocity  and  failure  diameter.  Good 
resistance  data  were  obtained;  however,  improper  voltage 
offset  on  the  recording  instrumentation  permitted  recovery  of 
only  the  later  portions  of  the  data. 

BE87-121.  This  experiment  was  a  repeat  of  BE87-118.  Voltage 
offset  was  corrected  and  good  resistance  data  were  obtained. 
In  loading  the  test  article,  however,  there  was  a  shortage  of 
the  20/80  mix  and  150  mm  of  30/70  mix  was  required  to  complete 
filling  of  the  driver  section.  Erratic  behavior  of  the 
detonation  velocity  in  the  conical  section  was  probably  caused 
by  the  long  section  of  30/70  mix  in  the  driver. 


Table  6.  Cylindrical  charge  test  matrix. 
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(a)  Plate  dent  test  at  41.28-mro  diameter  not  included 

(b)  Plate  dent  test  at  101.6-nun  diameter  not  included 


Table  7.  Conical  charge  test  matrix. 
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Table  8.  Big  Eagle  experiments. 
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BE87-240.  This  experiment  was  a  repeat  of  BE87-118  and  -121. 
In  this  experiment,  erratic  resistance  probe  data  were 
observed.  This  was  probably  caused  by  air  bubbles  that  were 
inadvertently  introduced  when  problems  were  encountered  during 
charge  loading  operations.  A  good  estimate  of  failure 
diameter  was  made,  however,  by  measuring  the  position  of  the 
mechanical  short  in  the  resistance  probe  which  was  recovered 
posttest . 

BE87-114.  No  pin  data  were  recovered.  No  explanation  could 
be  found  for  this  loss  of  data. 

In  most  cases,  charges  were  initiated  with  an  RP-80  and  8  mm  of 
Oeta-Sheet.  For  mixes  of  30/70  and  less  (less  HMX) ,  an  additional  layer 
of  several  millimeters  of  40B/60  was  added  for  additional  boosting. 

Good  ionization  pin  signals  were  observed  for  mixtures  of  40/60  and 
greater.  Adequate  signals  were  observed  for  30/70  mixes,  but  erratic 
signals  were  observed  in  the  initial  20/80  experiments.  In  later  20/80 
experiments,  shorting  pins  which  rely  on  mechanical  shorting  from 
pressure  effects*  rather  than  on  ionization  effects,  were  used  to 
register  arrival  time  of  the  detonation  front. 

DATA  ANALYSIS  METHODS 

A  brief  review  of  the  analysis  methods  used  to  reduce  the  data  will 
be  presented.  This  review  will  include  the  analysis  for  TOA,  resistance 
probe,  streak  camera  and  plate  dent  data. 

TQA  Analysis.  By  making  four  TOA  measurements,  six  velocities  can 
be  calculated.  The  detonation  velocity  results  in  Table  8  indicate  the 
average  of  the  six  velocities  with  the  maximiim  and  minimum  variations 
about  the  average  indicated  by  the  ±  values.  In  some  cases  only  three 
TOAs  were  recorded,  resulting  in  only  three  calculated  velocities.  As 
indicated  previously,  the  uncertainty  in  any  given  calculated  value  of 
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detonation  velocity  based  on  TOA  measurements  is  less  than  0.5%. 
Representative  ionization  pin  data  are  shown  in  Figure  7. 


Resistance  Probe  Analysis.  Voltage -time  data  were  recorded  on  the 
Lecroy  system.  The  change  in  position  of  the  shock- induced  mechanical 
short  in  the  resistance  probe.  Ax,  is  directly  related  to  the  change  in 
voltage,  AV,  through  the  relation 

Ax  -  AV/(kI)  (2) 

where  k  is  the  wire  resistance  per-unit- length  (0.2809  O/mm)  and  I  is 
the  current.  The  voltage- time  data  may  be  differentiated  at  various 
positions  to  determine  the  variation  of  detonation  velocity  with  x. 
Detonation  failure  is  indicated  by  an  abrupt  change  in  slope  of  the 
voltage -time  data  or  by  first  indication  of  constant  voltage.  Figure  8 
shows  the  results  of  a  resistance  probe  experiment,  BE87-118.  The  data 
were  analyzed  and  processed  using  the  Lecroy  CATALYST  program  on  an  IBM 
PC. 


Independent  TOA  measurements  are  made  at  several  positions  to 
establish  a  precise  position- time  reference  for  the  resistance  probe 
data.  Since  the  position  of  the  pins  and  the  times  at  which  the 
detonation  front  arrives  at  the  pins  are  precisely  known,  a  direct 
relation  can  be  established  between  the  recorded  voltage  and  pin 
positions.  Without  this  information,  AV/At  could  still  be  accurately 
determined  but  the  position  along  the  charge  axis  at  which  that  value  of 
AV/At  occurred  could  not  be  precisely  determined. 

From  the  geometry  of  the  conical  charges  the  diameter  at  any 
position  along  the  axis  is  precisely  known  and  is  calculated  from  the 
relation 


d  -  di  -  0.05237  x  (3) 

where  d|  is  the  charge  diameter  at  the  entrance  to  Che  conical  section 
and  X  is  the  distance  from  the  beginning  of  the  conical  section  to  the 
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37 


Figure  7.  Ionization  pin  data,  BE87'112. 
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Figure  8.  Resistance  probe  data,  BE87-118. 


location  where  the  diameter  is  d.  The  coefficient  is  based  on  a  conical 
angle  of  3  degrees. 

1 1 

Craig  has  shown  that  an  overshoot  in  the  propagation  of 
detonation  along  the  charge  axis  of  up  to  one  charge  diameter  can  occur 
at  detonation  failure.  This  occurs  because  the  detonation  at  any  given 
station  along  the  axis  of  a  conical  charge  is  always  slightly 
overdriven.  Thus,  conical  charge  data  tends  to  indicate  failure  at  a 
diameter  that  is  slightly  smaller  than  the  true  failure  diameter.  The 
lower  bound  on  failure  diameter  determined  in  the  conical  charge 
experiments  has  been  adjusted  to  reflect  this  overshoot. 

Streak  Caaera  Photography.  The  length  calibration  along  the  axis 
of  the  charges  in  the  streak  camera  experiments  was  established  by 
making  a  still  photo  of  the  test  article  and  using  the  known  length  of 
the  test  article  (or  a  millimeter  scale  included  in  the  photo)  as  the 
reference.  The  time  base  was  established  through  the  precisely  known 
geometry  of  the  film  plane  radius  relative  to  the  rotor  axis  and  the 
rotor  speed.  Thus,  the  film  writing  speed,  W,  is  determined  as 

W  -  4jrrn  (4) 

where  r  is  the  film  plane  radius  and  n  is  the  rotor  speed  in  rev/sec. 
Writing  speed  as  determined  above  is  accurate  to  within  0.25%.  The  film 
was  read  using  a  modified  fiche  reader.  The  effects  of  detonation 
overshoot  at  failure  are  considered  in  determining  the  failure  diameter 
and  equation  (3)  applies  as  before.  ^ 

Plate  Dent  Analysis.  A  simple  empirical  technique  has  been 
12 

developed  by  Smith  to  obtain  a  rough  determination  of  CJ  pressure  in 
an  explosive.  The  experimental  setup  consists  of  a  cylindrical  charge 
positioned  in  the  center  of,  and  normal  to,  a  steel  plate  which  has  been 
calibrated  against  known  explosives.  The  CJ  pressure  for  the  explosive 
is  determined  by  measuring  the  depth  of  the  dent  produced  by  detonation 
of  the  charge.  The  relationship  between  dent  depth  and  CJ  pressure  is 
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"cj  -  "  y 


(5) 


where  F  is  a  factor  which  is  related  to  the  charge  diameter,  y  is  the 
dent  depth  in  millimeters  and  is  in  kbar.  The  values  of  F  used  in 
the  two  plate  dent  experiments  (BE87-241  and  -242)  are  F  -  34.81  for  a 
diameter  of  41.28  mm  and  F  -  8.67  for  a  diameter  of  101.6  mm. 

BE87-241  was  conducted  on  the  40/60  mix  at  a  diameter  of  41.28  mm. 
A  dent  depth  of  6.73  mm  was  measured  which  results  in  a  calculated  CJ 
pressure  of  234  kbar  (23.4  GPa) .  This  is  in  good  agreement  with  the  CJ 
pressure  calculated  with  TIGER  of  237  kbar  (23.7  GPa).  The  good 
agreement  between  experiment  and  calculation  may  be  due  to  the  fact 
that,  at  this  diameter,  the  40/60  mix  is  approaching  ideal  behavior  and 
the  heavy  HMX  loading  makes  the  BKWR  parameters  appropriate  for  "RDX- 
like”  explosives. 

BE87-242  was  conducted  primarily  as  a  go,  no-go  experiment  for  the 
10/90  mix  at  a  diameter  of  101.6  mm.  A  witness  plate  was  used  to 
confirm  detonation  and  provide  an  estimate  of  CJ  pressure.  Dent  depth 
was  measured  in  this  experiment  at  8.59  mm  which  results  in  a  calculated 
CJ  pressure  of  74  kbar  (7.4  GPa).  This  is  well  below  the  theoretical  CJ 
pressure  of  180  kbar  (18.0  GPa)  calculated  with  TIGER.  This  suggests 
that  significant  diameter  effects  affected  the  results  of  this 
experiment  even  at  the  large  diameter.  In  addition,  the  10/90  mix  may 
be  more  "composite -like"  than  "RDX-like",  in  which  case  the  BKWR 
parameters  are  less  appropriate  than  for  the  40/60  case  above. 

DISCUSSKW 

Detonation  Velocity-DlaBeter-Hlxture  Analysis.  The  results  in 
Tables  8  and  9  are  plotted  in  Figures  9  through  12.  In  Figure  9 
detonation  velocity  as  a  function  of  inverse  charge  diameter  is  plotted 
for  the  five  HHX/emulsion  mixtures  for  which  data  were  obtained.  In 
Figure  10  the  logarithms  of  the  absolute  values  of  the  slopes  of  the 
linear  portions  of  the  lines  of  constant  mixture  in  Figure  9  are  plotted 
as  a  function  of  percentage  HMX  as  a  means  for  better  understanding  the 
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Figure  10.  Plot  of  logarithm  of  absolute  value  of  A  (slope  of 
D  vs  1/d  at  large  diameters  from  Figure  9)  vs  percentage  HMX. 


Figure  11.  Detonation  velocity  as  a  function  of  percentage 
for  diameters  of  12.55  nun,  26.10  mm  and  51.50  mm. 


PERCENTAfiE  HMX 

Figure  12.  Definition  of  detonation  failure  envelope. 
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trends  in  the  data.  In  Figure  11  detonation  velocity  as  a  function  of 
percentage  HMX  is  plotted  for  the  three  diameters  for  which  data  were 
obtained.  In  Figure  12  the  boundary  between  regions  of  detonation  and 
detonation  failure  is  estimated  by  plotting  the  results  of  the 
cylindrical  and  conical  charge  experiments  in  charge  diameter-percentage 
HMX  space. 

Figure  9  provides  a  comprehensive  description  of  the  results  of  the 
cylindrical  and  conical  charge  experiments.  Detonation  velocity  is 
plotted  as  a  function  of  inverse  charge  diameter  for  each  HMX/emulsion 
mixture.  The  open  symbols  are  data  from  Table  8  for  the  cylindrical 
charges.  The  solid  symbols  are  points  read  from  streak  camera  film 
records  and  resistance  probe  data  for  the  conical  charges.  The 
horizontal  bars  are  estimates  of  failure  diameter  based  on 
interpretation  of  the  streak  camera  film  records  and  resistance  probe 
data.  In  the  streak  camera  film  records,  indications  of  possible 
detonation  failure  are  a  sudden  change  in  slope  of  the  record  or  a 
sudden  change  of  light  intensity.  In  the  resistance  probe  data, 
indications  of  possible  detonation  failure  are  a  sudden  change  in  the 
slope  of  the  voltage -time  record  or  the  sudden  occurrence  of  constant 
voltage.  Uhile  the  results  in  some  cases  contain  large  uncertainties  in 
failure  diameter,  the  results  in  Figure  12  appear  somewhat  exaggerated 
because  they  are  plotted  as  inverse  charge  diameter. 

Figure  9  shows  that,  for  a  given  mixture,  a  straight  line  through 
the  cylindrical  charge  data,  extrapolated  to  the  infinite  diameter 
(TIGER)  velocity  represents  the  results  very  well.  This  is  consistent 
with  the  behavior  of  numerous  molecular  explosives  as  shown  by  Campbell 

13 

and  Engelke  It  is  reasonable  to  assume  that  for  diameters 
significantly  larger  than  failure  diameter,  there  is  a  consistancy,  or 
as  stated  by  Campbell  and  Engelke,  a  "regularity,"  in  the  relationship 
between  slopes  of  constant  mixture  and  percentage  HMX.  These  slopes  in 
Figure  9  were  measured  and  the  logarithm  of  their  absolute  values 
plotted  against  percentage  HMX  in  Figure  10.  Because  the  detonation 
velocities  of  the  40/60  mixture  appeared  to  be  anomalously  low  relative 
to  the  other  mixtures,  a  straight  line  fit  to  that  data  in  Figure  9  was 
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not  made.  However,  the  slope  of  a  line  through  that  data  was  measured 
and  those  results  included  in  Figure  10,  the  purpose  being  to  confirm 
the  suspected  "irregularity"  of  that  data. 

Some  discretion  has  been  exercised  in  deciding  how  best  to  fit  the 
data  in  Figure  10.  Giving  preference  to  the  20,  30  and  50  percent  HMX 
data,  a  curve  has  been  fit  by  hand,  recognizing  that  the  slope  of  the 
100  percent  HMX  mixture  is  very  nearly  zero  and  thus  influences  the 
behavior  of  the  curve  for  mixtures  greater  than  60  percent  in  the  manner 
shown.  For  mixtures  of  less  than  20  percent  HMX  the  curve  has  simply 
been  extrapolated  to  the  pure  emulsion  where  the  only  constraint  on  the 
result  is  that  the  slope  must  be  finite.  Results  in  this  region  can 
only  be  treated  as  an  estimate. 

The  results  in  Figure  10  clearly  show  that  the  40/60  mixture  data 
are  deviant.  However,  the  results  can  also  be  used  to  provide  an 
estimate  of  what  the  slope  of  the  40/60  mixture  should  be  for 
consistency.  This  estimate  is  indicated  in  Figure  9  with  a  dashed  line. 
In  addition,  the  results  of  Figure  10  have  also  been  used  to  provide  an 
estimate  of  the  slopes  for  10  percent  HMX  and  the  pure  emulsion.  These 
results  are  also  shown  in  Figure  9  as  dashed  lines. 

In  principle  it  is  possible  to  use  the  results  shown  in  Figure  9 
(horizontal  bars)  to  determine  the  detonation  velocity  at  failure  for 
each  mixture.  However,  due  to  the  large  negative  slopes  at  and  near 
failure,  and,  in  some  cases,  the  large  uncertainties  in  estimates  of 
failure  diameter,  it  is  not  possible  to  provide  a  definition  of  the 
failure  boundary  with  any  degree  of  confidence.  A  general  observation 
can  be  made,  however,  that  for  low  percentage  HMX  mixtures  failure 
occurs  at  large  diameters  at  low  velocities,  while  for  high  percentages 
HMX  mixtures  failure  occurs  at  small  diameters  at  high  velocities.  This 
result  is  in  qualitative  agreement  with  the  behavior  indicated  in  Figure 
1. 


The  short  lines  crossing  the  20/80  and  30/70  data  in  Figure  9 
indicate  that  at  those  diameters  detonation  was  not  observed  in 
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cylindrical  charge  experiments.  Since  the  results  of  the  20/80  and 
30/70  conical  charge  experiments  indicate  detonation  at  smaller 
diameters,  the  adequacy  of  the  boosters  used  to  initiate  these  charges 
becomes  suspect. 

Figure  11  shows  variation  of  detonation  velocity  with  percentage 
HMX  for  constant  diameter.  The  results  of  the  TIGER  calculations  (Table 
5)  are  included.  As  shown,  charges  with  higher  percentage  HMX  and 
larger  diameters  result  in  increased  detonation  velocity  as  well  as  a 
closer  approach  to  the  TIGER  results  which  are  assumed  to  be  a  good 
estimate  of  detonation  velocity  as  infinite  diameter.  The  open  symbols 
are  plots  of  the  results  from  Table  8  for  the  cylindrical  charges.  The 
solid  symbols  are  plots  of  points  determined  from  Figure  9  which  could 
be  determined  with  good  confidence.  Because  of  the  limited  data  and 
large  uncertainties  in  velocity  at  failure,  no  attempt  has  been  made  in 
Figure  11  to  establish  the  detonation  failure  boundary.  As  noted 
before,  qualitative  agreement  with  the  behavior  described  in  Figure  1  is 
demonstrated. 

In  Figure  12  an  estimate  of  the  detonation  failure  envelope  was 
made  by  constructing  a  curve  between  the  regions  of  detonation  and 
detonation  failure.  The  curve  shown  was  constructed  assuming  that  100 
percent  HMX  fails  at  less  than  1-mm  diameter  and  that  this  limit 
dominates  the  behavior  of  the  curve  for  high-percentage  HMX  mixtures. 
The  curve  is  only  approximately  defined  and  additional  experiments  would 
be  required  to  precisely  define  the  boundary,  especially  for  mixtures 
where  the  HMX  content  is  less  than  40  percent.  The  40/60  data  have  been 
included  here  even  though  prior  analysis  indicated  that  the  measured 
detonation  velocities  are  anomalously  low  relative  to  the  other 
mixtures .  Preference  has  been  given  to  the  conical  charge  data  over  the 
cylindrical  charge  data  in  constructing  the  curve.  Note  that  two 
cylindrical  charge  data  points  indicating  failure  lie  well  within  the 
detonation  region.  These  are  the  same  two  points  indentlfled  in  Figure 
9  which  may  have  been  inadequately  boosted.  Questions  about  adequacy  of 
initiation  may  apply  to  other  cylindrical  charge  data  as  well.  It  is 
felt  that  the  curve  in  Figure  12  is  reliable  to  within  ±  IS  percent  for 
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mixtures  containing  more  than  20  percent  HMX.  For  mixtures  containing 
less  than  20  percent  HMX  the  curve  can  be  created  only  as  an  estimate. 
A  precise  determination  of  the  failure  for  pure  emulsion  would  remove 
much  of  the  uncertainty  in  this  region. 

Diameter  Effect  Analysis.  The  diameter  effect  analysis  of  Campbell 

13 

and  Engelke  was  applied  to  the  50/50  mixture  data  in  Figure  9. 
Campbell  and  Engelke  determined  that  the  following  functional  form  was 
descriptive  of  the  behavior  of  a  large  number  of  heterogeneous 
explosives  in  a  cylindrical  geomety, 

D-D^{1-A/(R-R^)1  (6) 

where  is  Che  detonation  velocity  at  infinite  diameter,  R  is  the 
charge  radius  and  A  and  R^  are  length  parameters .  The  50/50  mixture  was 
selected  for  this  analysis  because  the  data  were  sufficient  to  provide  a 
reasonable  definition  of  the  downward  concavity  of  the  D  vs  1/d 
relationship  near  detonation  failure.  Four  points  on  the  D  vs  1/d  curve 
from  Figure  9  were  selected  and  fit  by  equation  (6).  The  results  are 
shown  in  Table  10  where  the  values  of  D/D^^  calculated  from  the  fit  are 
compared  with  the  experimentally  determined  values.  As  shown,  a 
reasonable  fit  is  achieved. 

If  equation  (6)  is  descriptive  of  all  other  HMX/emulsion  mixtures 
(there  is  no  reason  to  think  that  it  isn't)  and  the  Campbell  and  Engelke 
regularities  are  present  in  the  data  (they  should  be),  then  a 
generalized  description  or  model  of  the  diameter  effect  for  HMX/emulsion 
mixtures  should  be  achievable  by  defining  the  dependence  of  A  and  R^  on 
percentage  HMX.  This  generalized  model  will  not  be  pursued  here  but 
should  be  addressed  in  future  analysis. 

HMX  Particle  Size  Effects.  Some  remarks  are  in  order  concerning 
the  40/60  and  40B/60  mixtures.  In  the  revised  test  matrix  all  mixtures 
utilized  only  the  coarse  grade  HMX.  The  role  of  the  40B/60  mixture  in 
the  revised  test  matrix  was  to  demonstrate  (in  a  very  preliminary  way) 
the  effect  of  fine  grade  HMX  on  detonation  performance.  The  rationale 
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in  addressing  this  objective  became  somewhat  confused  when  the 
detonation  velocities  of  the  40/60  mixture  were  determined  to  be 
anomalously  low.  Undoubtedly,  failure  diameter  and  the  D  vs  1/d 
behavior  (Figure  9)  are  affected  and  the  meaningfulness  of  a  comparison 
with  the  40B/60  data  becomes  questionable.  As  indicated  in  Table  9, 
there  is  no  clear  difference  in  diameter  at  failure  between  the  two 
(LE87- 127/128  and  LE87-137),  though  the  actual  differences  may  be  lost 
in  the  uncertainties.  One  does  not  conclude  that  the  failure  diameters 
are  the  same,  nor  does  one  conclude  that  they  are  not. 

Also  note  that,  if  the  40B/60  data  from  Table  8  are  plotted  in 
Figures  9,  10  and  11  (this  was  not  done),  excellent  agreement  with  the 
other  data  occurs.  Based  on  the  Campbell 'Engelke  model  of  diameter 
effects,  it  would  probably  be  appropriate  to  include  those  data  in 
Figures  9,  10  and  11,  since  they  would  lie  on  the  nonconcave  portion  of 
the  curve  described  by  Equation  (6),  i.e.,  the  dashed  line  in  Figure  9. 
A  simple  examination  shows  that  these  data  do  lie  on  that  line. 

At  this  point,  we  terminate  discussion  of  the  40/60  and  40B/60  data 
in  that  nothing  further  can  be  contributed  to  real  understanding  of  the 
detonation  performance  of  HHX/emulsion  mixtures. 


i 
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sEcntm  5 
OMICLDSIONS 

At  the  completion  of  the  experimental  program  described  in  this 
report,  it  was  clear  that  many  additional  experiments  would  be  required 
to  provide  a  complete  description  of  the  detonation  performance  of 
HMX/emulsion  mixtures.  Given  that  it  was  not  possible  to  obtain  good 
mixes  for  HMX  loadings  greater  than  50  percent,  the  additional 
experiments  needed  are  for  detonation  velocity-diameter  data  for 
mixtures  below  20  percent  HMX  and  for  more  refined  data  on  failure 
diameter  and  velocity  at  failure  for  all  mixtures.  Schedule  and  funding 
constraints  did  not  permit  additional  experiments.  Nevertheless,  it  was 
possible  to  provide  velocity-mixture-diameter  data  over  a  region 
sufficient  to  provide  a  definitive  description  of  the  behavior  of 
HMX/emulsion  mixtures  over  the  range  of  20  to  SO  percent  HMX  and  a 
qualitative  description  over  the  entire  range  of  mixtures.  The  results 
of  this  research  should  be  valuable  in  pointing  the  way  for  future 
research  on  the  detonation  performance  of  molecular/emulsion  explosive 
systems . 

The  main  conclusions  of  this  research  are  summarized  as  follows; 

(1)  At  a  given  diameter,  an  increase  in  HMX  loading  results  in  an 
increase  in  detonation  velocity.  As  HMX  loading  increases,  the  failure 
diameter  decreases  and  the  detonation  velocity  at  failure  increases. 
The  trend  in  velocity  at  failure  was  evident  in  the  data  but  precise 
measurements  of  velocity  at  failure  were  not  obtained.  The  qualitative 
behavior  predicted  in  Figure  1  has  been  confirmed. 

(2)  The  diameter  effect  on  detonation  velocity  for  HMX/emulsion 
mixtures  is  adequately  described  by  the  Campbell -Engelke  model.  This 
was  clearly  demonstrated  for  the  50/50  mixture  and  the  trend  of  downward 
concavity  on  the  D  vs  1/d  curve  was  observed  for  other  mixtures. 

(3)  Addition  of  approximately  20  percent  HMX  to  the  emulsion 
results  In  a  substantial  Increase  In  Initiation  sensitivity.  This 
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observation  is  based  on  the  premise  (which  is  believed  to  be  true  for 
the  HMX/emulsion  mixtures)  of  an  inverse  relationship  between  failure 
diameter  and  initiation  sensitivity.  An  interesting  area  of  future 
research  would  be  to  investigate  the  combined  effects  of  particle  size 
and  mass  loading  of  a  molecular  explosive  on  initiation  sensitivity  and 
failure  diameter. 

(4)  The  particular  emulsion  formulation  developed  in  this  research 
exhibited  (in  most  cases)  remarkable  physical  stability.  After 
preparation  of  the  emulsion,  HMX  was  added  in  a  high  shear  mixing 
operation  at  room  temperature,  and  the  mixture  was  subjected  to 
additional  shear  during  preload  and  test  article  loading  operations  - 
also  at  room  temperature.  Samples  of  the  pure  emulsion  have  been 
observed  for  over  one  year  with  only  slight  clouding  of  the  emulsion. 

(5)  A  variety  of  physical  properties  of  the  HMX/emulsion  mixtures 
were  observed,  depending  on  the  HMX  loading,  HMX  particle  size  and  shear 
history.  The  four  classes  of  properties  observed  were: 

(a)  Stiff  grease  -  obtained  for  HMX  loadings  from  10  to  50 
percent  for  coarse  grade  HMX. 

(b)  Noncohesive  powder  •  obtained  for  mixtures  of  coarse  grade 
HMX  for  loadings  above  SO  percent  and  for  fine  grade  HMX  at  a  loading  of 
50  percent. 

(c)  Hard  •  obtained  for  coarse  grade  HMX  at  a  loading  of  50 
percent  when  subjected  to  high  shear. 

(d)  Nontacky  putty  -  obtained  for  the  40B/60  mixture. 

The  castable  (hard)  and  putty- like  consistencies  have  properties  that 
make  them  of  Interest  In  a  number  of  potential  applications.  Because 
the  physical  nature  of  the  50/50  cast  mixture  Is  unknown.  Its  Initiation 
sensitivity  to  shock,  friction  and  spark  effects  should  be  thoroughly 
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characterized  during  the  initial  phase  of  any  future  research  program 
involving  that  material. 

(6)  While  no  firm  conclusions  are  made  based  on  the  40B/60  and 
40/60  data  relative  to  the  effects  of  HMX  particle  size  on  detonation 
performance,  the  downward  concavity  observed  in  the  D  versus  1/d  curve 
demonstrates  that  the  hot-spot  mechanism  controls  propagation  of 
detonation  at  small  diameters;  and  thus  it  would  be  predicted  that  a 
decrease  in  particle  size  would  result  in  a  decrease  in  failure  diameter 
as  well,  since  for  a  given  mass  loading  the  number  of  initiation  sites 
would  increase  in  direct  proportion  to  the  inverse  cube  of  the  particle 
diameter. 
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10  '  PROGRAM  IS  "BALANCE".  CALCULATES  OXYGEN  BALANCE 
20  '  FOR  EXPLOSIVES  COMPOSED  OF  CARBON,  HYDROGEN,  OXYGEN, 

30  '  NITROGEN,  SODIUM  AND  ALUMINUM.  CALCULATIONS  ARE  BASED 
40  '  ON  REDUCTION  OF  ALL  REACTANTS  TO  H20  AND  C02. 

50  ' 

60  '  WRITTEN  BY  J.  RENICK,  AIR  FORCE  WEAPONS  LABORATORY,  24  FEB  88. 
70  ' 

80  DIM  C(10) ,H(10) ,0(10) ,AL(10) ,N(10) ,NA(10) ,OB(10) .X(10) 

90  DIM  A$(10) 

100  MWC-12 
110  MWH-1 
120  MWN-14 
130  MWO-16 
140  MWAL-27 
150  MWNA-23 

160  INPUT  "EXPLOSIVE  NAME?";Z$ 

170  LPRINT  Z$ 

180  LPRINT  “  " 

190  LPRINT  "  " 

200  INPUT  "NUMBER  OF  REACTANTS  -";NR 

210  FOR  I-l  TO  NR  STEP  1 

220  INPUT  "REACTANT  NAME?":A$(I) 

230  INPUT  "CARBON-" ;C(I) 

240  INPUT  "HYDROGEN-" :H( I) 

250  INPUT  "OXYGEN-"; 0(1) 

260  INPUT  "ALUMINUM-" ;AL(I) 

270  INPUT  "SODIUM-" :NA( I) 

280  INPUT  "NITROGEN-" :N( I) 

290  MWT-C(I)*MWC+H(I) *MWH+N ( I ) *MWN+0 ( I ) *MWO+ AL ( I ) *MWAL+NA ( I ) *MWNA 

300  OB(I)— (2*C(I)+H(I)/2+1.5*AL(I).0(l))*(MW0/MWT) 

310  NEXT  I 
320  LPRINT  "  " 

330  LPRINT  "  " 

340  LPRINT  "REACTANT/C/H/N/O/AL/NA/OXYGEN  BALANCE" 

350  LPRINT  "  " 

360  FOR  I-l  TO  NR  STEP  1 

370  LPRINT  A$(I) ;C(I) :H(I) ;N(I) ;0(I) :AL(I) :NA(I) :OB(I) 

380  NEXT  I 

390  FOR  I-l  TO  NR  STEP  1 
400  PRINT  "  " 

410  PRINT  A$(I) 

420  INPUT  "NASS  FRACTION  -";X(I) 

430  NEXT  I 
440  OXB-0 

450  FOR  I-l  TO  NR  STEP  1 
460  OXB-OXB  +  X(I)*OB(I) 

470  NEXT  I 
480  LPRINT  "  " 
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